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Oxygen adatoms at SrTiOs(OOl): A density-functional theory study 
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We present a density-functional theory study addressing the energetics and electronic structure 
properties of isolated oxygen adatoms at the SrTiO3(001) surface. Together with a surface lattice 
oxygen atom, the adsorbate is found to form a peroxide-type molecular species. This gives rise to a 
non-trivial topology of the potential energy surface for lateral adatom motion, with the most stable 
adsorption site not corresponding to the one expected from a continuation of the perovskite lattice. 
With computed modest diffusion barriers below 1 eV, it is rather the overall too weak binding at both 
regular SrTiO3(001) terminations that could be a critical factor for oxide film growth applications. 
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I. INTRODUCTION 



Apart from its use in photo-catalytic and sensing 
applications^, the SrTiO3(001) surface is also receiving 
increasing attention as a suitable substrate material for 
thin film growth 3 -^. For the latter context the numerous 
reported surface reconstructions^' 7 -'^ ' 10 ' 11 ' 12 , partly in 
sensitive dependence of the applied annealing tempera- 
ture, indicate a complex surface kinetics, which needs 
to be understood and controlled when aiming at growth 
experiments tailored to the atomic-scale. This view of 
the surface being far from equilibrium under the em- 
ployed experimental condition a 12 ' 13 is corroborated by 
first-principles thermodynamic calculation a 13 ' 14 ' 15 , which 
identified the unreconstructed surface as the only equi- 
librium termination. In the resulting focus on the ad- 
sorption kinetics, a first important step is to establish 
detailed insight into the nature of the bond and the con- 
comitant binding and diffusion properties of all surface 
species involved. Studying the interaction of oxygen with 
SrTiO3(001) is in this respect a natural starting point, in 
light of the important role played by oxygen exposure for 
surface preparation and growth or gas-sensing applica- 
tions. Seemingly simple at first glance, preceding experi- 
mental and theoretical work has already revealed intrigu- 
ing aspects even of the initial surface intermediates. Not 
untypical for wide band-gap oxides 3 -^, spectroscopic and 
kinetic investigations suggest the existence of a molecu- 
lar peroxide "Og 2 " species at the surfac e - 17 ' 18 . This is 
confirmed by recent density-functional theory (DFT) cal- 
culations by Piskunov et al^, which demonstrate that 
this peroxide species results at both regular SrTiO3(001) 
terminations actually from the bonding of an O adatom 
to a surface O ion of the oxide lattice. 

The formation of such a dioxygen moiety may readily 
give rise to a non-trivial topology of the potential en- 
ergy surface (PES) for lateral diffusion of adsorbed O 
atoms across the surface, with direct consequences for 
the role of O adatoms in the growth kinetics. Aiming 
to elucidate the latter we therefore employ DFT to ex- 
plicitly map this PES for isolated O adatoms at both 
regular SrTiO3(001) terminations and discuss the impli- 
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FIG. 1: (Color online) Schematic illustration of the crystal 
structure of cubic SrTiOa and its two regular (001) termi- 
nations, featuring a Ti02 or SrO composition in the topmost 
layer. O atoms are depicted as small, light gray (blue) spheres, 
Ti atoms as medium-size, gray (green) spheres and Sr atoms 
as large, dark (red) spheres. 



cations of the obtained energetic data with respect to the 
presence and mobility of this surface species. The com- 
puted surface electronic structure indeed suggests that 
the complex formed between the O adatom and lattice 
oxide ion may be characterized as a molecular peroxide 
species. While this surface bond is at both terminations 
strong enough to shift the most stable adsorption site 
away from the one that would correspond to a continua- 
tion of the perovskite lattice, the absolute bond strength 
is only about the same as in gas-phase O2 and the dif- 
fusion barrier with 0.8 eV quite modest. This points to 
the overall stabilization, rather than the mobility of a 
corresponding O adatom species as critical factor for the 
growth kinetics at this surface. 



II. THEORY 



Identified as 

-13.14,15 



the only equilibrium surface 
structures^ii^i'ii, the regular Ti0 2 - and SrO- 
terminations illustrated in Fig. [1] form a natural basis 
for our fundamental study addressing the energetics and 
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electronic structure of O adatoms at cubic SrTiO3(001). 
All corresponding DFT calculations were performed 
with the CASTEP 20 code using a plane-wave basis 
together with ultrasoft pseudopotentials 2 ^ as provided 
in the default library, and the GGA-PBE functional 22 
to treat electronic exchange and correlation (xc). Tests 
using the LDA did not yield any qualitative differences 
with respect to the conclusions put forward below. 
The surfaces were modeled with supercell geometries, 
using inversion-symmetric slabs based on the GGA-PBE 
optimized bulk lattice constant (ao = 3.938 A), with O 
adsorption at both sides and with a vacuum separation 
exceeding 13 A. Apart from the central three layers of 
the slab, all substrate atomic positions were fully relaxed 
until the absoulute value of all ionic forces dropped 
below O.leV/A, with further tightening of this criterion 
not yielding any significant changes with respect to 
surface geometry and concomitant energetics. This 
relaxation was equally applied to the z-height of the O 
adatom, with its x and y coordinates frozen at different 
positions to map the PES for lateral adatom motion. 

The centrally targeted energetic quantity is the binding 
energy of an O adsorbate, denned as 

Eb = ^ [-^OQsurf — -Esurf — ^0 2 (gas)] , (1) 

where £b<a SU rf is the total energy of the O covered sur- 
face, -Esurf the total energy of the clean surface, £ , o 2 (gas) 
the total energy of the gas-phase O2 molecule (all three 
computed at the same plane- wave cutoff), and the fac- 
tor 1/2 accounts for the fact that adsorption is at both 
sides of the slab. The binding energy is thus referenced 
to a free O2 molecule in its spin-triplet ground state (cal- 
culated in a (18 x 18 x 18) A supercell), and a negative 
value of £b indicates that adsorption is exothermic. Both 
-Eb@surf and E sur f were furthermore consistently com- 
puted in a non spin-polarized way, with systematic tests 
confirming that spin-polarization at the surface does not 
play any significant role for the description of all relevant 
low-energy parts of the PES. 

Systematic tests with one O adatom at different high- 
symmetry sites in a (1 x 1) surface unit-cell verified that 
i?b is numerically converged with respect to the basis 
set to within ±20meV when employing a plane wave ki- 
netic energy cutoff of 430 eV. Due to long-range geomet- 
ric relaxations 2 ^ achieving a similar convergence with re- 
spect to the number of surface unit cells in the slab model 
is more challenging. This concerns prominently those ad- 
sorption sites where the O adatom protrudes significantly 
from the surface, e.g. on-top of the lattice Ti-atom on 
the Ti02 termination. However, at the same time such 
sites are energetically so unfavorable that they do not 
even belong to the saddle point regions of the adatom 
PES that are relevant for the diffusion barriers. Con- 
sequently, we employed (3 x 3) surface unit-cells with a 
(2x2x1) Monkhorst-Pack grid 2 ^ for the reciprocal space 
integrations. Test calculations in larger surface unit-cells 
indicate that this set-up provides a very good description 



of the isolated adatom limit with respect to the most sta- 
ble site and diffusion barriers addressed here. The slab 
models comprised in total seven and 11 atomic layers for 
the Ti02- and SrO-termination, respectively, which en- 
sure that all energetically relevant parts of the adatom 
PES discussed below are converged to at least within 
±50meV. 



III. RESULTS 
A. Geometric structure and energetics 

Figure [2] summarizes the computed PES for lateral O 
adatom motion at both regular terminations. As sus- 
pected, the shape is in both cases non-trivial, with the 
global minimum, i.e. most stable adsorption site, not 
corresponding to the one expected from a mere continu- 
ation of the perovskite lattice. The geometric adsorption 
structure at this most stable site shown in Fig. [3J for 
both terminations is indeed very consistent with the an- 
ticipated reason^ for this PES topology in form of strong 
bonding between the adsorbate and a SrTiO3(001) sur- 
face oxygen anion (vide infra). At the SrO-termination 
this directly coordinated lattice oxygen atom is pulled 
out of the surface while the underlying Ti atom, denoted 
as Ti(l) in Fig. [3J is pressed into the bulk material, 
thereby stretching the Ti-0 bond from 1.97 A to 2.20 A. 

On the other hand the distance to the adsorbed O a d 
atom is quite short, and with 1.52 A rather character- 
istic for a molecular peroxide "Oj" 2 " species&i 6 .. With 
respect to the surface normal the dioxygen complex is 
tilted with a Ti(l)-0-O a( j angle of 120° and points be- 
tween two surface Sr atoms. However, the barrier for a 
90° rotation of the complex around the z-axis to point 
between an adjacent pair of Sr atoms is only 0.24 eV, cf. 
Fig. O One could thus view the dioxygen moiety as a 
tilted, singly-bonded 77 1 "end-on" configuration with re- 
spect to the underlying Ti(l) cation in the second layer, 
with the tilt direction exhibiting some dynamics at finite 
temperature. 

At the Ti02-termination only atoms in the topmost 
layer are significantly affected by the presence of the oxy- 
gen adatom. The oxygen-oxygen bond measures 1.47A 
and is therewith again in the ballpark expected for a 
molecular peroxide species. It is tilted towards one of 
the adjacent Ti surface atoms denoted as Ti(l) in Fig. 
[31 such that the distance between O a d and the latter is 
1.91 A and the Ti(l)-0-O ad tilt angle is 62°. This geom- 
etry is thus more reminiscent of an asymmetric, doubly- 
bonded i] 2 "side-on" configuration with respect to the 
coordinated Ti(l) cation, which is in fact the more usual 
form in dioxygen-Ti complexes 2 ^ and with the respec- 
tive bond-lengths and angles of the Ti(l)-0-O a d moiety 
agreeing again specifically with those of peroxotitanium 
complexes 26 . Despite this coordination a flipping motion 
of the dioxygen moiety to the adjacent Ti atom denoted 
as Ti(2) in Fig. [3J has only a small computed barrier of 
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FIG. 2: (Color online) Potential energy surface PES(a;,j/) for lateral motion of an O adatom at the SrO- and Ti02-termination, 
upper and lower panels respectively. The lateral coordinates x and y are given in units of the bulk lattice constant ao = 3.938 A, 
with the shown range for each termination exactly corresponding to the surface unit-cell definition shown in Fig. [T] At the 
SrO-termination, the position (0,0) is thus directly atop a surface Sr atom, while at the Ti02-termination the position (1/2, 
1/2) is atop a surface Ti atom. The white diamonds denote the actually calculated points within the irreducible wedge and 
neighboring isolines correspond to a difference of 0.25 eV. 



0.30 eV as shown in Fig. [21 which again indicates some 
dangling dynamics at finite temperature. 

Even in this most stable site at the Ti02-termination 
the computed binding energy is with 0.30 eV still en- 
dothcrmic with respect to gas-phase O2. However, due 
to the known shortcomings of the employed gradient- 
corrected xc functional such absolute binding energy val- 
ues have to be considered with care. In this respect, it 
is quite gratifying that the preceding study by Piskunov 
and coworkers^, which employed a hybrid functional in- 
cluding some non-local Hartree-Fock exchange, reports a 
rather similar binding energy of 0.66 eV at the same most 
stable site and with nearly identical geometry parame- 
ters. In fact, when redoing our calculations in the smaller 



(2 x 2) surface unit-cell employed by them we arrive at 
a virtually coinciding £4, = 0.60 eV. This seems to indi- 
cate that a proper account of the long-range geometric 
relaxations enabled in the larger surface unit-cells accessi- 
ble with the computationally much less demanding semi- 
local functional have a larger effect on the binding ener- 
getics than the improved xc description introduced by the 
hybrid functional. Furthermore, we believe that the ad- 
ditional GGA advantage of being able to perform a much 
more thorough exploration of the PES has direct conse- 
quences for the comparison of the results obtained for the 
SrO-termination. Within the understanding of the PES 
topology shown in Fig. [2] it appears as if the restricted 
testing of three possible high-symmetry adsorption sites 
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FIG. 3: (Color online) Perspective view of the most stable O 
adsorption geometry at the SrO- and TiCh-termination, up- 
per and lower panel respectively. Atoms specifically referred 
to in the text are labeled. O atoms are depicted as small, light 
gray (blue) spheres, Ti atoms as medium-size, gray (green) 
spheres and Sr atoms as large, dark (red) spheres. 



has misled the authors of ref*^ to erroneously identify the 
transition-state for the rotational motion of the oxygen 
moiety around the z-axis as most stable adsorption site, 
while the PES area around the above described most sta- 
ble site of lower symmetry was not scanned in their study. 
Recomputing this transition state in the smaller surface 
unit-cell we again achieve excellent agreement between 
our GGA-PBE and their reported hybrid functional en- 
ergetics, which fortifies our confidence in the here cal- 
culated slightly exothermic E\> = — 0.20eV at the most 
stable site shown in Fig. [3l 

At least at the Ti02-termination the binding is thus 
clearly too weak to generate O adatoms through disso- 
ciative adsorption at the ideal terraces. This is com- 
monly expected for insulating oxides 3 and has been ver- 
ified explicitly for related materials like SnC>2 27 , MgO 28 
and TiOz^. Considering the high vacancy formation 
energies^, it would, however, energetically clearly be 



possible to create an O adatom by adsorbing an oxy- 
gen molecule into a surface O vacancy at both regular 
terminations. Once established at the surface, the com- 
puted barriers for diffusive motion to the adjacent bind- 
ing site of 0.81 eV and 0.68 eV at the SrO- and Ti0 2 - 
termination respectively, cf. Fig. [51 demonstrate that 
the 0-O a d bond is then much more readily broken com- 
pared to the desorption of the 0-O a( j molecule so that O 
adatoms should be a relatively mobile species. 



B. Electronic structure 

The preceding section has discussed the geometric 
structure and adsorption energetics within the picture of 
a quasi-molecular species formed between the O adatom 
and a directly coordinated lattice O anion. This view is 
nicely supported by the adsorbate-induced change of the 
electron density shown in Fig. 0J This quantity, com- 
monly referred to as "difference electron density"—, is 
obtained by subtracting from the electron density of the 
O@SrTiO 3 (001) system both the electron density of the 
clean SrTiO3(001) surface and that of an isolated oxygen 
atom. Here, the atomic positions of SrTiO3(001) and of 
the O atom are taken to be same as those of the relaxed 
adsorbate system. In this way, the presentation high- 
lights the charge rearrangement upon bond formation, 
which as clear from Fig. U primarily takes place between 
the O adatom and its lattice O partner. 

Some further quantification of this adsorbate-induced 
charge redistribution can be obtained from an inspection 
of the Mulliker* 3 ^ and Hirshfeld 32 charges on the vari- 
ous surface species. With the two charge partitioning 
schemes yielding the usual differences in absolute num- 
bers, they nevertheless unanimously predict the same 
trend in that the charge surplus at the lattice O anion 
is essentially evenly distributed over the formed dioxy- 
gen moiety after adsorption. In other words, while the 
charges on all other atoms in the topmost two layers are 
virtually unaffected by the adsorbed O atom, the charge 
on the directly coordinated lattice O anion gets about 
halved with the O adatom exhibiting a charge of roughly 
equal magnitude. For the Mulliken (Hirshfeld) scheme 
this means specifically that at the Ti02-termination the 
computed charge of — 0.73e (— 0.33e) at the surface O 
anion, diminishes to — 0.49e (— 0.16e) upon adsorption 
with the O adatom exhibiting a similar charge of — 0.38e 
(— 0.16e). At the SrO-termination the equivalent num- 
bers are -0.83e (-0.38e) and -0.54e (-0.20e) for the 
surface anion before and after O adsorption, as well as 
— 0.58e (— 0.30e) for the O adatom. The thereby sug- 
gested formation of a molecular peroxide "0^~ 2 " species 
is corroborated by an analysis of the density of states 
(DOS) at the surface, which we focus for the present 
purposes on the O 2s dominated states around 20 eV be- 
low the bulk valence band maximum. Figure O compiles 
the corresponding calculations and compares them to the 
oxygen-projected DOS of bulk Sr02, a material in which 
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FIG. 4: (Color online) Difference electron density An(r) plot highlighting the formation of a dioxygen moiety at the SrO- and 
TiCh-termination, left and right panel respectively. Shown is the charge rearrangement within the plane perpendicular to the 
surface defined by the O a( j adatom and its directly coordinated lattice O and Ti(l) partners, cf. Fig. [3] with red (blue) areas 
with solid (dashed) contour lines depicting charge depletion and accumulation, respectively. Neighboring isolines indicate a 
difference of 10~ 3 e/A . The superimposed spheres indicate the position of atoms within the plane, using the color and size 
scheme for the individual species as in Figs. [l]and|3] 



0^ 2 -type bonding occurs naturally. Apparently, O ad- 
sorption at SrTiC>3(001) leads at both terminations to 
two sharp fingerprint states above and below the regular 
SrTiC>3 O 2s group, which almost perfectly match the 
energetic position of the equivalent states in bulk SrC>2. 
This confirms the interpretation in terms of of a molecu- 
lar peroxide species, the bond of which leads to the non- 
trivial PES topology for O adatom diffusion described in 
Section 3. A. 



IV. CONCLUSIONS 

In conclusion we have presented a density-functional 
theory study analyzing the potential-energy surface for 
lateral diffusion of O adatoms at both regular termina- 
tions of SrTiO3(001). The O adatom is found to form 
a quasi- molecular "0^" 2 "-type species with a directly co- 
ordinated surface O anion of the substrate lattice. This 
bond is strong enough to shift the most stable adsorption 
site away from the one expected from a mere continuation 
of the perovskite lattice, and gives rise to a non-trivial 
PES topology that allows for some dangling dynamics 
of the dioxygen moiety at finite temperature. In addi- 
tion, the computed barriers for diffusive hops to adjacent 
binding sites are at both terminations well below 1 eV, so 
that mobility of O adatoms is not expected to represent 
a bottleneck at typical film growth conditions. On the 



other hand, the overall bond strength is at both termina- 
tions only about equal to the one of a O2 molecule, and 
specifically slightly exothermic at the SrO-termination 
and slightly endothermic at the TiCvtermination with 
respect to this gas-phase reference. At least at the latter 
termination the binding is thus clearly too weak to gen- 
erate O adatoms through dissociative adsorption at ideal 
SrTiO3(001) terraces. In view of the high O surface va- 
cancy formation energies, an at least energetically much 
more plausible alternative would be that O adatoms are 
created by adsorbing an oxygen molecule into such a de- 
fect. Once established at the surface, the computed mod- 
est diffusion barriers show that the "0^~ 2 " peroxide-type 
bond is broken rather readily, suggesting that such O 
adatoms would provide a mobile species that could an- 
neal further defects or attach to created islands. 
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FIG. 5: (Color online) Density of states (DOS) analysis sup- 
porting the interpretation of the formed O a d-0 complex in 
terms of a molecular peroxide species. The topmost panels 
show the total DOS of the topmost two surface layers of the 
SrO-termination of clean SrTiO3(001) before (a) and after 
(b) O adsorption. The next two panels show the DOS pro- 
jected onto the 2s (green line) and 2p (blue line) states of 
the O a d atom and its coordinated lattice O partner for the 
SrO- (c) and the Ti02-termination (d), identifying in par- 
ticular two sharp states above and below the O 2s group as 
fingerprints for the formed O a( j-0 moiety. Panel (e) shows 
the same projected DOS of bulk Sr02 exhibiting essentially 
the same characteristics. The bulk valence band maximum 
(as approximately reached in the central slab layer) is taken 
as zero reference in all panels. 
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